function f=IsPowerPositive(A)
for n=1:10
    AP=A^n;
    if (min(min(AP))>0)
        f=true;
        return;
    end;
end
f=false;
return;